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Abstract 

Iterated Bernstein polynomial approximations of degree n for continu- 
ous function which also use the values of the function at i/n, i = 0,1, ... ,n, 
are proposed. The rate of convergence of the classic Bernstein polynomial 
approximations is significantly improved by the iterated Bernstein poly- 
nomial approximations without increasing the degree of the polynomials. 
The close form expression of the limiting iterated Bernstein polynomial 
approximation of degree n when the number of the iterations approaches 
infinity is obtained. The same idea applies to the g-Bernstein polynomials 
and the Szasz-Mirakyan approximation. The application to numerical inte- 
gral approximations which gives surprisingly good results is also discussed. 
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1 Introduction 

The Bernstein polynomials [l| have been used for approximations of functions 
in many areas of mathematics and other fields such as smoothing in statistics 
and constructing Bezier curves [see for examples] which have important 

applications in computer graphics. One of the advantages of the Bernstein poly- 
nomial approximation of a continuous function / is that it approximates / on 
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[0, 1] uniformly using only the values of f at i/n, i = 0,1, ... ,n. In case when the 
evaluation of / is difficult and expensive, the Bernstein polynomial approximation 
is preferred. 

The properties of the Bernstein polynomial approximation have been studied 
extensively by many authors for decades. However the slow optimal rate 0{l/n) 
of convergence of the classical Bernstein polynomial approximation makes it not 
so attractive. Many authors have made tremendous efforts to improve the per- 
formance of the classical Bernstein polynomial approximation. Among many 
others, Butzeryi introduces linear combinations of the Bernstein polynomials 



and Phillips [5] proposes the g-Bernstein polynomials which is a generalization of 
the classical Bernstein polynomial approximation. However, Butzerj^'s approx- 
imation involves not only the Bernstein polynomials of degree n but also degree 
of 2n which requires more sampled values of the function to be approximated at 
the 2n + 1 rather than n + 1 uniform partition points of [0, 1]. The g-Bernstein 
polynomial approximates a function / only when g ^ 1. For g > 1, it seems that 
f{z) has to be an analytic complex function on disk {z : < r}, r > g, so 
that the g-Bernstein polynomial approximation of degree n has a better rate of 
convergence, 0{q~'^), than the best rate of convergence, of the classical 

Bernstein polynomial approximation of degree n [seeQ, Q, for example]. If g > 1, 
the g-Bernstein polynomial approximation of degree n uses the sampled values of 
the function at n + 1 nonuniform partition points of [0, 1]. These points except 
t = 1 are attracted toward t = when g is getting larger so that the approxi- 
mation becomes worse in the the neighborhood of the right end-point. This is a 
serious drawback of the g-Bernstein polynomial approximation which limits the 
scope of its applications. 

In this paper, we propose a simple procedure to generalize and improve the 
classical Bernstein polynomial approximation by repeatedly approximating the 
errors using the Bernstein polynomial approximations. This method involves only 
the iterates of the Bernstein operator applied on the base Bernstein polynomials 
of degree n and the sampled values of the function being approximated at the 
same set of ?t, + 1 uniform partition points of [0, 1]. The improvement made by the 
g-Bernstein polynomial approximation with properly chosen g can be achieved by 
the iterated Bernstein polynomials without messing up the right boundary. 
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2 Preliminary Results About the Classical Bern- 
stein Polynomial 

Let / be a function on [0, 1]. The classical Bernstein polynomial of degree n is 
defined as 

n 

Mj(t) = M^^^f{t) = J2f{i)Bm{t), o^t<:i, (1) 

i=0 

where B„ is called the Bernstein operator and Bni(t) = i = 0, . . . ,n, 

are called the Bernstein basis polynomials. Note that the Bernstein polynomial 
of degree n, Bn ■*/, uses only the sampled values of f at tni = i/n, i = 0,1, ... ,n. 
Note also that for i = 0, . . . ,n, 

f3^^{t) = {n + l)Bm{t), O^t^l, 

is the density function of beta distribution beta{i + l,n + 1 — i). Let Yn{t) be a 
binomial b{n,t) random variable. Then E{Yn{t)} = nt, var{y„(t)} = E{Yn{t) — 
nty = nt{l-t), E{Yn{t)-ntY = nt{l-t){l-2t), and MJ{t) = E[f {Ynit) /n}]. 
The error of f is 

Err{l«/}(t)=B«/(t)-/(t). (2) 

Let / be a member of C^''^[0,1], the set of all continuous functions that have 
continuous first r derivatives. C[0, 1] = C^°''[0, 1]. Let the modulus of continuity 
of the rth derivative /'•^•' be 

u,{5)= max _ ^ > q. 

\s-t\<5 

About the rate of convergence of B^^^/ we have the following well known results 
[see 8. 

Theorem 1. Suppose f E C^'^''[0, 1], r = 0, 1. For each n > 1 

|Err{B«/}(t)| = |B„/(t) - ^ CrU-^/'urin-'/'), 

where Cr is a constant depending on r only. One can choose Cq = 5/4 and 
Ci = 3/4. 

n 

The result according to r = is due to Popoviciu[9i]. The order of approxima- 
tion of / e C^^^ [0, 1] by arbitrary polynomials is given by the theorem of Dunham 



Jackson 10 
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Theorem 2 (Dunham Jackson). Suppose f G C^^''[0, 1], r ^ 0. For each n > r 
there exists a polynomial Pn of degree at most n so that 

|P„(t)-/(t)| 

where C'^ is a constant depending on r only. Ifr = 0, one can choose Cq = 3. 

11| about the asymptotic formula 



The following is a result of Voronovskaya 
of the Bernstein polynomial approximation. 

Theorem 3 (E. Voronovskaya). Suppose that f has second derivative f". Then 



Err{BiiV}W=B„/(t)-/(t) 



\^nt) + U^{t\ 
In n 



(3) 



where Snit) is a sequence of functions which converge to as n 



oo. 



From Theorem [3] it follows that the best rate of convergence of bI'V, 

as 

n — *• oo, is 0{n~^) even if / has continuous second or higher order derivatives j^. 
This is not as good as in the case of arbitrary polynomial approximation in which 
if / has continuous rth derivative then the rate of convergence of a sequence of ar- 
bitrary polynomials P„ of degree at most n can be at least o{n~^) 10|]. Bernstein 



12l | generalizes this asymptotic formula to contain terms up to the (2/c)th deriva- 



tive and proposes a polynomial constructed based on both f{i/n) and f"{i/n), 
i = 0,1, . . . ,n. Butzer 0] considers some combinations of Bernstein polynomials 
of different degrees and shows that they have better rate of convergence which is 
much faster than 0{l/n). Costabile et al 13| generalize the linear combinations 



of the Bernstein polynomials proposed by of jj] , ij] and 15| . The g-Bernstein 



poly„om,als of fl has better rate of eor;verge„ce HoweveT if < , < 1. the 

g-Bernstein polynomials of function / do not approximate /. For g > 1, the 
g-Bernstein polynomials do approximate / at a rate of C(g~'") but f{z) has to 
be analytic in a complex disk with radius greater than g. The analyticity of / 
may be too restrictive for applications. Even if we are sure that / is analytic, 
we have to deal with the choice of g. In some cases, the approximations are very 
sensitive to the choice of g. 



3 The Iterated Bernstein Polynomials and the 
Rate of Convergence 

The error Err{]Bi^'*/}(t) is also a continuous function on [0,1] whose values at 
ti = i/n, i = 0,1, ... ,n, depend on f{ti), i = 0,1, ... ,n, only. So we can 
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approximate this error function by the Bernstein polynomial Bi"^^[Err{Bi^''/}](t) 
and then subtract the approximated error function from Mn^ f{t) to obtain the 
second order Bernstein polynomial of degree n 



Kt^fit) = Bi^VW - B^^^[Err{Bi,^7}](t). (4) 

This idea is closely related to, although was not initiated by, the proposal of 
Bernstein 1^ in which the second derivative rather than the error of the Bernstein 
polynomial is approximated. Inductively, 

B(f+i)/(t)=B(f)/(t)-B„,{B(f)/(t)-/(t)}, k^l. (5) 

This iteration procedure can be performed further until a satisfactory approxi- 
mation precision is achieved because the error Err{Bn*'''/(t)} = M^'^ f{t) — f{t) 
can be estimated by B„{Bi^VW - /W) = - B^f+^VW- 

Lemma 4. Generally the k-th order Bernstein polynomial of degree n can he 
written as 



/(t), O^t^l. (6) 



Define B°/(t) = f{t). Then the error of the k-th Bernstein polynomial of degree 
n can be written as 

Err{B(f)/(t)} = BW/(t) -/(t) = Q (-l)^-^/^ = -(I-B„)V(t), (7) 



where I = B° is the identity operator. 
Proof. 

B('=+^)/(t) = BW/(t) - B„{B(f)/(t) - fit)} 



irX^7W+B„/(t) 



i=l ^ ^ 1=1 ^ / 

g Q i-iy-'MUit) + E (, ! i) i-^y-'Km + 
fe+i 



T.(''^,}(-^y"w)- (8) 



i=l 



By induction, ([T]) and ([8]) assure that ([6]) is true for every positive integer k. 
Equation ([7]) is then obvious. □ 
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The limit of B^/(t), as — > oo, has been given by Kehsky and Rivhn 16 
A short and elementary proof of [la|'s result is given by [l3]. After we finished 
the first version of this paper, we realized that 18| obtained the formula ([7]) 
and investigated the properties of B^f''/(t) using simulation method. The cost of 
B^f''/(t) is only some simple algebraic calculations in addition to the evaluation 
of f at i/n, i = 0,1, ... ,n. 

About the iterates of the Bernstein operator we have the following result. 

Lemma 5. For k ^ 1, 

n 

Kfit)=Y.f&n-\B^m. k^l- O^t^l, (9) 
i=0 

where B°i?„j(t) = Bniif), and 

B^+i5„,(t) = B4B^5„,}(t), k^l. (10) 

When k = 1, 

n 

MlBUt)=Y.Mi)Bnj{t). (11) 

j=0 

Proof. The theorem can be easily proved by induction and the fact that the 
Bernstein operator is linear. □ 

By ([6]) and ([7]) we have 

Theorem 6. The k-th Bernstein polynomial approximation can he calculated 
inductively as 

" /h.\ 

4'=V(t) = 5^/(^)$^ . (-l)^-X-^5™(t), k^l, O^t^l. (12) 

i=0 i=l ^-^^ 

Clearly, for every k ^ 1, Mn^ preserves linear functions. Therefore 

Err{BW/(t)} = {/(^) - /W} E H ("ir 'Br^5„(t), k^l;0^t^l. 

1=0 j=l ^"'^ 

(13) 

Expression (fT2i) can easily implemented in computer languages using iterative 
algorithm. Define indicator functions 
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Then B^i{t) = B„J„(t) = Bi'^4,(t), i = 0,l,...,n, and, by Theorem El 
and ( fT3l) can be simphfied as 



Mi^^f{t) = J2f{iP^n^Ut), k^l, O^t^l. (15) 

i=0 
n 

Err{B(f)/(t)} = J] {/(^) - /(t)} B(f)j„,(t), k^l;0^t^l. (16) 

w = K'^ut) = E (^) (-i)^-^Bri5„,(t). (17) 

i=i ^-^^ 

The following theorem shows that the iterated Bernstein polynomials , like the 
classical ones, have no error at the endpoints of [0, 1]. 

Theorem 7. For any function f defined on [0, 1] and any integer k 0, 

4'V(0) = /(O), B(f)/(1) = /(1). (18) 

Proof. It is known that B°5„i(t) = 5„i(t) = /„i(t) for t = 0, 1, z = 0, . . . , ra. 
Assume that B^~^_B„j()f:) = /„i(t) for t = 0, 1, i = 0, . . . , n, and some k 1. By 
Theorem [9l if t = 0, 1, 

n n 
M'^Br,i{t) = J2Bm{iPn~\Br,,m = B„,, (i) /„, (t) = 5„,(t) = 4,(t). 

i=o i=o 

So by induction, for all nonnegative integers k, t = 0,1, and i = 0,...,n, 
M'^B^iit) = Ini{t). By (HID, we have, if t = 0, 1, 

j=i Vj/ 

Thus by if t = 0,1, 

n 

^^n^f{t) = Y.f{^Imit)=f{t). 



i=Q 



□ 



Clearly, for each k ^ 1, B^fV(i) 

can be written as 
B(f)/(t) = )S„(t) 
where i^i*^'* = {f^^)ix(n+i) is an (n + 1) row vector, and 

Bnit) = {Bnoit), . . . , Bnnit)V- 
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If = 1, 

Define (n + 1) x (n + 1) square matrix *8n = Bn{u'^) = (&ij)(n+i)x(n+i) where 

"""" ^ (n' n' • • • ' n) ' 

That is 

bij = Bn,i-i{^-^), i,j = l,...,n + l. 

It is easy to see that 23„ is nonsingular and have all the eigenvalues in (0, 1] 
among them exactly two are ones which correspond to eigenvectors tt„ and l„+i = 
(1, . . . , 1)^ G -R""*"^. We have the following theorem. 

Theorem 8. For k ^ 1, 

1=1 ^ ^ 

where = In+i, the {n + l)st order unit matrix. If k ^ 1, 

pik^i) ^ Fit\ln+i - + Fi^). (20) 

Proof. It is easy to show that 

p(2) _ 2^(1) _ I7^(i)m 

By induction, (fT9|) and (!20|) can be easily proved. □ 

More importantly, we have 
Theorem 9. The "optimal" Bernstein polynomial approximation of degree n is 
^^rT^fit) = Ft^B^S) = F«!8„^B„(t), (21) 

where 

Fi-) = hm = FW!8;\ (22) 

fe^oo 

Moreover, mI^^ preserves linear functions. 

Proof. Since all the eigenvalues of matrix *B„ are in (0, 1] and exactly two of 
them are ones, all the eigenvalues of matrix In+i — are in [0, 1) and exactly 
two of them are zeros. Thus limfc^oo(-^n+i — ^n)^ = O, the zero matrix. Because 
Mn^ preserves linear functions for any positive integer k, so does This can 

also be proved by the following facts that 

fW!8„ = F« if and only if F^^^^ = F« 



and that Fi^^^„ = F^^^ is true provided that / is hnear. 



□ 
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Numerical examples (see ^ show that the maximum absolute approximation 
error seems to be minimized by "optimal" Bernstein polynomial approximation 
Bn°°''/(t) if / is infinitely differentiable. For nonsmooth functions such as f{t) = 
\t — 0.5| and fixed n, it seems that the maximum absolute approximation error is 
minimized by the iterated Bernstein polynomial approximation Mn^ f{t) for some 
k. 

The next theorem shows that if A; > 1 then Mn^ f{t) is indeed a better poly- 
nomial approximation of / than the classical Bernstein polynomial. 

Theorem 10. Suppose that f G Cd^.^[0, 1], d/^r = 2(/c — 1) +r and r = 0, 1. Then 

|Err{B('=)/(t)}| = |B('=)/(t) - f{t)\ ^ Cln-^uj,Jn-'/'), (23) 
where C'l^ is a constant depending on r and k only. 

Proof. This result follows easily from Theorems [1] and [31 □ 

Remark 3.1. From this theorem with k = 2 and r = 0, we see that if f has 
continuous second derivative then the rate of convergence of the second Bernstein 
polynomial approximation Mn^ f is at least o{n~^). 

Remark 3.2. From Theorem [721 with k = 2 we see that if f has continuous 
fourth derivative, then the rate of convergence ofMn^f can be as fast as 0{n~'^). 
This seems the fastest rate that f can reach even if f has continuous fifth or 
higher derivatives. 

Remark 3.3. It can also be proved that if f has continuous {2k)th derivative, 
then the rate of convergence ofM^^f can be as fast as 0{n^^). Although these 
improvements upon Mnf{t) are still not as good as those stated in Theorem\^ 
they are good enough for application in computer graphics and statistics. 

Remark 3.4. It is a very interesting project to investigate the relationship be- 
tween C'l^ and k, and the rate of convergence ofM^^^'f which is conjectured to be 
exponential. 

4 The Derivatives and Integrals of Bi^^/(t) and 
Applications 

4.1 The Derivatives of B?V(^) 
Theorem 11. For any positive integers k and r ^ n, 

Jr I n~r k / ]c\ 
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where is the rth forward difference operator with increment h = 1/n, A/(t) 
f{t + h)-f{t), 



i=0 

Proof, li k = 1, it is well known that for any function / 

n-l 

„ M^f(t]=n 
dt ' ' dt 



= ^M^fit)=nJ2^f{i)Bn-iAt)- (25) 



i=0 



Assume that (HM with r = 1 is true for the kth iterated Bernstein polynomial of 
any function /. By we have 

|B(f+^)/(t) = I [B(f)/(t) - MM'^fit) - fit)}] 

= + |B„/(t) - |BjB(f)}/(t). (26) 

It follows from ([25]) and that 

, n— 1 

-BaBW}/(t) = ABi'=)/(^)i?„_i,(t) 

n— 1 n ^ / 1 \ 

= «EE/(i) E (* (-l)'-AB^'B„,(-)B„-„(() 

i=0 j=0 l=\ ^ ^ 

n— 1 fe / ; \ 

= -EE(j(-l)'"^®n/(^)5n-M(t). (27) 
1=0 ^=1 ^ ^ 

Combining ([2S]), ([2]), (El]), and ([27]) we arrive at 

|Br^)/(t) = n5^5^(-l)-^r + MABrVa)i?.-M(t). (28) 

i=0 7=1 \ J / 



The proof of (IMj) with r = 1 and /c ^ 1 is complete by induction. Similarly ([2i 
with r ^ 1 and k ^ 1 can be proved using induction. □ 



It is not hard to prove by adopting the method of [8| that 

Theorem 12. (i) If f has continuous rth derivative /*■''•' on [0, 1], then for each 
fixed k, as n ^ oo, ^Mn^ f{t) converge to f^'^\t) uniformly on [0, 1]. 
(a) If f in bounded on [0, 1] and its rth derivative f^'^\t) exists att E [0, 1], then 
for each fixed k, as n oo, -^Mn^ f{t) converge to f^^\t). 
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Numerical examples show that the larger the r is, the slower the above con- 
vergence is. 

For any positive integers k, the second derivative of the iterated Bernstein 
polynomial M^''^ f is 

|^B('=)/(t) = n{n , ^'(®rV) {Ti)Bn-2M- (29) 

4—n ,—1 / 



i=0 j=l 

T rn 11 fllOTI , 



It is well known that if / is convex on [0, 1], then ^B>^n^ f{t) ^ and thus Bi^V(i) 
is also convex and Mn^f{t) ^ f{t) on [0,1]. So the classical Bernstein polyno- 
mials preserve the convexity of the original function and has nonnegative errors. 
However examples of §H] show that when 2 the iterated Bernstein polynomial 
^i^) f does not preserve the convexity of the original function unconditionally. 
The iterated Bernstein polynomials still preserve the monotonicity of / if it is 
not too "flat" anywhere. 

Theorem 13. // / is strictly increasing (decreasing) on [0,1], for any k ^ 1, 
Mn^f{t) is also strictly increasing (decreasing) on [0,1]. 

Proof. The theorem is true for = 1 even if / is increasing (decreasing), but not 
strictly, on [0, 1]. It suffices to prove the theorem when / is strictly increasing 
on [0, 1]. Assume that the theorem is true for some k ^ 1. Since / is strictly 
increasing on [0, 1], B^/(t) are also strictly increasing on [0, 1] for all A; ^ 1. □ 

Remark 4.1. If k = 1, the condition of strict monotonicity is not necessary. 
However, if k > 1, the condition of strict monotonicity can be relaxed. For 
example, f{x) = x, if 0^x<l/3, = 1/3, if 1/3 ^ x < 2/3, and = x — 1/3, 
if 2/3 ^ X ^ 1. It can be shown that /(x) < for x in a neighborhood of 

X = 1/2. 

4.2 The Integrals of B?V(^) 

The following theorem is very useful for implementing the iterative algorithm in 
computer languages. 

Theorem 14. Suppose f is continuous on [0, 1]. For 1 ^ k ^ oo and x G [0, 1], 

we have 

/ BW/(t)cit = ffS^x) = Fi'^^S^x), (30) 
where 5„(x) = {Sno{x) , . . . , Snn{x)y and 

px px 
Sni{x) = / Bni(t)dt = / Pni{t)dt, S'„j(l) = 

Jo Jo 
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Corollary 15. If g is continuous on [a,b], a < b, then for 1 ^ k ^ oo, 

j\{t)dt - ^E/i? = (31) 

where Fn''^ is calculated based on f{t) = j;^g[ci + {b — a)t]. 

Remark 4.2. Note that numerical integration l[JJ\) does not involve any integrals. 
It contains only algebraic calculations. See ExamplelB of ^ for some numerical 
examples. 

Proof. The theorem follows immediately from Theorems [S] and M D 

The following theorem follows immediately from Theorems [10] and [T4l 
Theorem 16. Under the condition of Theorem UR for any x G [0, 1] 

^Cln-'-^uj,,M~'^% (32) 



K^fm- / fit)dt 
Jo 

where C'l^ is a constant depending on r and k only. 

5 Iterated Szasz Approximation and Iterated q- 
Bernstein Polynomial 

The idea used to construct the iterated Bernstein polynomial approximation is 
simple and very effective. The same idea seems also applicable to other operators 



or approximations such as the Szasz operator [19|] [or the Szasz-Mirakyan (Mi- 
rakja) operator] and the g-Bernstein polynomial with g > 1. We will give some 
numerical examples in ^ and the analogues of results of Section [3] could be be 
obtained by using the analogue results about the rate of convergence of the Szasz- 



Mirakyan approximation 20||. We hope these would inspire more investigations 
with rigorous mathematics. 

5.1 Iterated Szasz Approximation 

The so-called Szasz-Mirakyan approximation is defined as 

oo 

^nf{x) = J2fii)Pn^{x), XG[0,Oo), (33) 
i=0 
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where / is defined on [0, oo) and Pni{x) = e~^^{nxy/i\. Note that, for a; > 0, 
Pni{x) is the probabihty that Vn{x) = i where Vn{x) is the Poisson random vari- 
able with mean nx. Since the binomial probability Bni{t) can be approximated 
by Pni{t) for large n, the Szasz-Mirakyan approximation can be viewed as an 
extension of the Bernstein polynomial approximation. The error of §„/ as an 
approximation of / is 

oo 

Err(§„/)(a;) = SJ{x) - f{x) = /(^)Pn.(x) - f{x), x G [0,oo). (34) 

i=0 

Applying the Szasz-Mirakyan operator to Err(Sn/)(x), we have 

oo 

S„{Err(§„/)}(a;) = ^lf{x)-§J{x) = ^ /(^)§„P„,(a;) - §„/(x), x e [0,oo). 

(35) 

So we can define the second Szasz-Mirakyan approximation as 

§(fV(a:) =Sn/(a;)-§„{Err(§„/)}(x), a;G[0,oo). (36) 

Theorem 17. 

oo ^ /] \ 

Si^V(a:)=5^/(^)5^(. (-l)^-V„-^Pn.(x), xG[0,oo). (37) 

1=0 j=l ^-^^ 

Clearly, for every k ^ 1, preserves linear functions and therefore 

oo ^ / 1 \ 

Err{§«/(x)} = {/(^) - fix)} H (-l)^"^Sr^P™(t), k ^ 1; x E [0, oo). 

i=0 j=l ^ 

(38) 

Figure H] gives an example of the iterated Szasz approximations. 



5.2 Iterated g-Bernstein Polynomial 

Let a; be a real number. For any q > 0, define the g-number 

^ ifg^l; 



[X\q 



1-q 

X, if g = 1. 



If X is integer, then [x]q is called a g- integer. For g 7^ 1, the g-binomial coefficient 
(Gaussian binomial) is defined by 



1, r = 0; 

(l^qr)^l-qr-l)...(l-q) , i ^ ^ 

0, r > n. 



13 



So 



n 



n 



^ r ^ n, q > 0, 



where empty product is defined to be 1. Thus the ordinary binomial coefficient 
(") is the special case when q = 1. G. M. Phillips jB| introduced the g-Bernstein 
polynomial of order n for any continuous function f{t) on the interval [0, 1] 

n 

Qn = Y,f (t!"^) Qm{t), n = 1, 2, . . . , 

i=0 

tf'^ = ^, Qm{t)=('')ffi{l-tq^-'), . = 0,1, 



where 



n. 



Clearly, B„/(t) = Q„i/(t) which is the classical Bernstein polynomial of order 
n. It has been proved that if < g < 1 then Q„g/(t) does not approximate / 
and that if g > 1 and f{z) is analytic complex function on disk {z : |z| < r}, 
r > q, then Qnqfit) has better rate of convergence, than the best rate of 

convergence, 0{n^^), of ]B„/(t) [see 

aa 

for example]. 

Note that if g > 1 then points t-'''' = are no longer uniform partition 

points of the interval [0, 1]. For fixed n, limg^oo^l^^ = 0, i < n. So all ti except 
tn^ = 1 are attracted toward as g getting large. However, interestingly, the 
larger the g is in a certain range, the closer the g-Bernstein polynomial approxi- 
mation Qnqfit) to fit). For a given n, if g is too large, the g-Bernstein polynomial 
approximation Qnqf{t) becomes worse in the neighborhood of the right end-point. 

Similarly we have the iterated g-Bernstein polynomials 



t) 



1=0 j=l ^-^^ 

See Figure [5] for an example of the iterated g-Bernstein polynomials. Comparing 
Figures [T] and O we see that increasing g from 1 to 1.1 does improve the approxi- 
mation on [0, 1] except at points in the neighborhood of the right end-point. The 
approximation near the right end-point could be worsen by applying the iter- 
ated g-Bernstein polynomials. The improvement can be achieved by the iterated 
Bernstein polynomials without messing up the right boundary. 



6 Numerical Examples 

In this section some numerical examples are given with the hope of more inves- 
tigations on the proposed methods with rigorous mathematics. 
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Example 1. Figure U\ shows the first three iterated Bernstein polynomials of 
f{t) = sin(27rt) and the errors where n = 30. The "optimal" Bernstein polynomial 
approximation is also plotted which seems to have almost no error. 

Example 2. Figure IB shows the first three iterated Bernstein polynomials of 
f{t) = sign(t — 0.5) (t — 0.5)2 differentiate hut not twice differentiahle function) 
and the errors where n = 30. The "optimal" Bernstein polynomial approximation 
is not plotted which becomes very had near the two endpoints. 

Example 3. Figure shows the first three iterated Bernstein polynomials of 
f{t) = |t — 0.5| and the errors wheren = 30. The "optimal" Bernstein polynomial 
approximation is not plotted which becomes very bad near the two endpoints. 

Example 4. Figure\^shows the first three iterated Szasz approximation of f{x) = 
0.25xe~^/2^ X ^ 0, and the errors where n = 10. 

Example 5. Figure shows the first three iterated q-Bernstein polynomials of 
f{x) = sin(7rx) and the errors where n = 30, g = 1.1. The performance of the 
approximation near t = 1 is very sensitive to q. 

Example 6. FigurelB shows the first three iterated Bernstein polynomials of the 
following function f{t) = \t — 0.5| and their derivatives where n = 30. 

Example 7. Figure shows the first three iterated Bernstein polynomials of the 
following function fit) and their derivatives where n = 30, 



This a convex function which has continuous first derivative but does not have a 
continuous second derivative. 



where /o(t) = v — \/r'^ — {t — uY is portion of a circle with radius r (a larger 
positive number) and centered at {u,v), u,v > 0, Pk(t) is a polynomial of degree 





k = 3. 



k 



Pk(t) — a^if — a^kt^ + a^^k-it'' ^ + ■ ■ ■ + cikit + Ofco- 



i=0 
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Table 1: Some results of numerical integrals {n = 5) 





k 






1 


5 


oo 


Exact value 


TT sm{'Kx)dx 


1.611471 


2.005416 


1.999203 


2 


lo e^dx 


1.746528 


1.718369 


1.718282 


1.718282 


Jo fix)dx 


0.3371903 


0.3413510 


0.3413443 


0.3413447 



Table 2: Some results of numerical integrals (n = 10) 





k 






1 


5 


oo 


Exact value 


TT sm{7rx)dx 


1.803203 


2.000146 


2.000000 


2 


lo e^^dx 


1.732389 


1.718285 


1.718282 


1.718282 


Jo f{x)dx 


0.3392624 


0.341345 


0.3413447 


0.3413447 



// we choose 

_ -30ts ± ^900tj - 40(25t^ - 



u 



Vr 



20 

then /(O) = /o(0) = 0, f{ts) = foits) = -3t5- We also have 



t — u „„, , 



Choose the coefficients of so that /(I) = Yl'i=o^ki = and the jth (j = 
0,1, . . . , k — 1) derivative at ts satisfy 

f^'\U) = j^jj^.a,,tr^ = fil\t\-'). 
i=j 

If r is large enough, say r = 70, 6 = 0.05, then f{t) is strictly convex and has 

(2) 

continuous positive second derivative f", but Mn f is still not convex because its 
second derivative is negative at some points near t = 0.4 (see Figurel^. 

Example 9. In the following Tables U\ andl^ we summarize some the results of 
numerical integrals on [0, 1] using our proposed method given in Corollary [731 for 
functions f{x) = 7rsin(7ra;), f{x) = , and f{x) = (p{x) = (l/-\/27r) exp(— x^/2). 



From these examples and the figures we see that the error is reduced sig- 
nificantly by using the iterated Bernstein polynomial approximation without in- 
creasing the degree of the polynomial. For non-smooth function, the maximum 
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error is reduced more than 50% by the third Bernstein polynomiaL It is also 
seen from Figure [3] that unhke the classical Bernstein polynomial approximation 
the iterated Bernstein polynomial approximation M^!^^ f seems not to preserve the 
convexity of / for > 1 in this case when / is not smooth. So it is necessary for 
Mn^f to preserve the convexity of / that / is smooth and /" is not too close to 
zero. For applications in numerical integrals and computer graphics, sometimes it 
is even much more expensive to evaluate the function / than the simple algebraic 
calculations. So it is significant to apply the iterated or the "optimal", if / is 
infinitely differentiable, Bernstein polynomial approximation. 
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Figure 1: The iterated Bernstein polynomials and errors when f{t) = sin(27rt). 
The error is minimized by Bi°°^/. 
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Figure 2: The iterated Bernstein polynomials and errors when /(t) = sign(t — 
0.5)(t — 0.5)^ which is differentiable on [0, 1] but not twice differentiable at t = 0.5. 
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Figure 3: The iterated Bernstein polynomials and errors when f{t) = |t — 0.5 1 
which is not different iable at t = 0.5. 
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Figure 4: The iterated Szasz approximations and errors when /(x) = 0.25xe~^/^, 
X ^ with n = 10. 



22 




Figure 5: The iterated g-Bernstein polynomials and errors when f(t) = sin(27rt) 
with n = 30, g = 1.1. 
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Figure 6: The iterated Bernstein polynomials and their derivatives when f{t) = 
|t — 0.5| which is convex but not different iable at t = 0.5. 
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Figure 7: The iterated Bernstein polynomials of f{t) as in Example [7] and their 
derivatives where f{t) is convex, differentiable on [0, 1] but not twice differentiable 
at t = 0.5. 
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Figure 8: The iterated Bernstein polynomials of / as in Example M and their 
derivatives. The function / is strictly convex but bI^V is not convex. 
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